The Lyman Continuum Polarization Rise in the QSO PG 1222+228Q 



Gregory A. Shields 
Department of Astronomy, University of Texas, Austin, TX 78712 

ABSTRACT 

o" ' , 

Some QSOs show an abrupt, strong rise in polarization at rest wavelength ~ 750 A. 
If this arises in the atmosphere of an accretion disk around a supermassive black hole, 
it may have diagnostic value. In PG 1222+228, the polarization rise occurs at the 
wavelength of a sharp drop in flux. We examine and reject interpretations of this 
feature involving a high velocity outflow. The observations agree with a model involving 
several intervening Lyman limit systems, two of which happen to coincide with the 
Lyman continuum polarization rise. After correction for the Lyman limit absorption, 
the continuum shortward of A912 is consistent with a typical power-law slope, a ~ —1.8. 
^vq This violates the apparent pattern for the Lyman limit polarization rises to occur only 

in "candidate Lyman edge QSOs". The corrected, polarized flux rises strongly at the 
wavelength of the polarization rise, resembling the case of PG 1630+377. The rise in 
■ polarized flux places especially stringent requirements on models. 
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1. INTRODUCTION 

Spectropolarimetric observations with the Hubble Space Telescope (HST) have revealed an 
unexpected rise in linear polarization in several QSOs (see review by Koratkar & Blaes 1999). 
These are radio quiet "candidate Lyman edge QSOs", in which the continuum flux drops rather 
rapidly at rest wavelengths A < 1000 A (Antonucci, Kinney, & Ford 1989; Koratkar, Kinney & 
Bohlin 1992). Models of accretion disk atmospheres predicted a reduced polarization in the Lyman 
continuum because of a diminished contribution of electron scattering to the opacity (Laor, Netzer, 
&; Piran 1990). With this motivation, Impey et al. (1995) and Koratkar et al. (1995) used the 
Faint Object Spectrograph (FOS) on HST to obtain ultraviolet spectropolarimetry of several QSOs 
with redshifts sufficient to bring the Lyman continuum within the observed wavelength band. The 
surprising result, in several cases, was a rapid rise in polarization in the Lyman continuum. From 
values <^ 1 percent in the optical and near ultraviolet, the observed polarization rises around rest 
wavelength 750 A to values ~ 5 percent in several objects, and to ~ 20 percent in PG 1630+377. 
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This phenomenon has inspired several attemps at explanation. Blaes and Agol (1996) found 
that, for effective temperatures T e ff ss 25, 000 K and low effective gravities, a polarization rise of up 
to ~ 5 percent at about the observed wavelength could occur naturally in QSO disk atmospheres. 
This results from the interplay of electron scattering, bound-free opacity, and the temperature 
gradient in the atmosphere. However, Shields, Wobus, and Husfeld (1998, hereinafter SWH) showed 
that the effects of the relativistic transfer function destroy the agreement between this model and 
observation. Beloborodov and Poutanen (1999) suggested a model involving Compton scattering 
in a corona or wind, but this model appears to have trouble giving the rapid rise in polarized flux 
observed in PG 1630+377 (Blaes & Shields 1999). Lee and Blandford (1997) discussed the possible 
role of scattering by resonance lines of heavy elements (see Section 5). 

SWH showed that, if the polarization is assumed to rise sharply at A912 in the rest frame of 
the orbiting gas, then relativistic effects would naturally produce the wavelength dependence of 
the observed polarization. This may offer a way of measuring the black hole spin, but the physical 
mechanism for the polarization rise remains unknown. 

PG 1222+228 is a B « 15.5 radio quiet QSO (Schmidt & Green 1983) whose polarization rise 
at A s=s 750 A coincides with a sharp drop in flux (Fig. 1, 2). Impey et al. (1995) noted this and 
attributed it to a coincidental Lyman limit system (LLS), corresponding to an identified absorption 
line system at z = 1.486. However, the coincidence of a broad absorption feature with a polarization 
rise also is observed for broad absorption line (BAL) QSOs. In these objects, outflowing gas at 
velocities ~ 10 4 km s _1 produces blueshifted absorption troughs, typically seen in the resonance 
lines of H I, C IV, N V, O VI, Si IV, and sometimes Mg II (Weymann et al. 1991; Arav, Shlosman, 
& Weymann 1997). Spectropolarimetric observations (e.g., Ogle 1997; Schmidt & Hines 1999; Ogle 
et al. 1999) often show a rise in polarization in the troughs, reaching values as high as ~ 8 to 10 
percent from ~ 1 percent at unabsorbed wavelengths. This is explained in terms of scattering of 
some of the continuum by an extended region that is not covered by the BAL flow (Hines & Wills 
1995; Goodrich & Miller 1995; Cohen et al. 1995). This pattern resembles the polarization rise 
and flux drop in PG 1222+228. 

This paper addresses two questions: (1) Does the polarization rise in PG 1222+228 result from 
an intrinsic absorber, analagous to the situation in the BAL QSOs? (2) If the drop in flux in PG 
1222+228 is an intervening LLS, what are the consequences of correcting the observed, polarized 
continuum for this absorption? 

2. INTRINSIC ABSORPTION IN PG 1222+228? 

We first consider the possibility that the flux drop and coincidental polarization rise in PG 
1222+228 results from some kind of intrinsic absorption. Two possibilities, considered below, are 
that it is a BAL outflow, or that it is an unusual, intrinsic LLS. 

In either case, one issue is the behavior of the polarized flux as a function of wavelength. As 
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discussed above, if the polarization rise results from the selective absorption of the directly viewed 
continuum but not the scattered continuum, one might expect a smaller drop (but generally not 
a rise) in the polarized flux, I p = pl\. In order to examine this, we have rebinned the data of 
Impey et al. (1995), kindly made available in reduced form by C. Impey and C. Petri (1999). 
These data consist of a spectrum with the G190H grating covering AA1575 to 2320 at 0.37 A per 
pixel, and a spectrum with the G270H grating covering the range 2224 to 3295 at 0.52 A per pixel. 
The G190H data shortward of A1994 have low signal-to-noise and were not presented by Impey 
et al. (1995). We used seven wavelength bins (in A): (1)1994-2224, (2) 2224-2287, (3) 2287-2319, 
(4) 2319-2492, (5) 2492-2761, (6) 2761-3029, (7) 3029-3295. Bins 2 and 3 involve an average of 
the two overlapping spectra, and bin 3 is a narrow bin containing the flux drop at A2300. The 
resulting values of I\ , q = Q/I and u = U/I are tabulated in Table 1, along with the polarization, 
p = \Jq 2 + u 2 , and its position angle, 9. For the polarized flux, we use the rotated Stokes flux Q' 
and polarization q' = Q'/I (cf. Koratkar etal 1995), referred to a position angle of 168 degrees. This 
is based on the mean polarization position angle of the HST data, which is in reasonable agreement 
with optical observations (Stockman et al. 1984; Webb et al. 1993). (The use of q' is appropriate 
if the position angle of the polarization is constant with wavelength. Table 1 supports this and 
also shows that the polarization p is reasonably consistent with q' .) These quantities are plotted in 
Figure 2. The shortest wavelength bin has larger polarized flux than the longer wavelength bins. At 
face value, this would weigh against an intrinsic absorber model for PG 1222+228; but it involves a 
single wavelength bin with substantial error bars. Therefore, we consider other aspects of the two 
outflow models. 



2.1. BAL Absorption 

The rest wavelength of the onset of the absorption feature in PG 1222+228 is ~ 750 A. Dome 
BAL QSOs show absorption by Ne VIII A775 (e.g., Arav et al. 1999; Telfer et al. 1999). This 
might be a candidate for the feature in question, in as much as BALs often set in at a wavelength 
somewhat blueshifted from the emission-line redshift. However, Ne VIII normally is accompanied 
by absorption in O VI A1035, N V A1240, and C IV A1550. There is no indication of broad C IV or 
Mg II absorption in the spectrum of PG 1222+228 (Sargent, Steidel, & Boksenberg 1988; Steidel & 
Sargent 1992). The HST spectrum shows a shallow trough at A3000 to A3070 that could be a weak 
O VI feature; but this may simply be a cluster of lines, including several strong La lines indentified 
by Impey et al. (1996). Photoionization models by Hamann (1997) indicate a range of ionization 
parameters for which the fractional abundance of Ne +7 exceeds that of +5 . However, given the 
normal ratio of oxygen to neon abundances, the O VI feature would likely be strong in a situation 
giving strong Ne VIII. 

The flux drop at A750 does not recover, with decreasing wavelength, in a way suggestive of a 
BAL (Figure 1). The spectrum has not fully recovered by rest wavelength 650 A, corresponding 
to an an outflow velocity of more than 40,000 km s _1 if attributed to Ne VIII. Some moderately 
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narrow "mini-BAL" features have been observed at such high velocities (Hamann et al. 1997), 
but true BAL troughs rarely reach such velocities. The same can be said in connection with the 
possibility that the A750 feature corresponds to a blend of features including N III, N IV, O IV, S 
VI, and Ne VIII seen in BAL QSO spectra at this wavelength (e.g., Arav et al. 1999). Moreover, 
given the range of ionization stages contributing to this blend, C IV and Si IV absorption would 
likely accompany it. 

Recent work has shown that BAL QSOs systematically have weak soft X-ray emission. BAL 
QSOs have optical to X-ray slopes a ox <; —2.0, whereas nonBAL QSOs tend to have a ox in the 
range -1.3 to -1.8 (Brandt, Laor, and Wills 1999). Here, a ox is defined by F X /F Q = (v x /v ) a ° x , 
where F x and F Q are the flux densities (F u ) at 2 keV and 3000 A, respectively. ROSAT pointing 
observations give a flux of 6 x 10~ 14 erg s _1 cm~ 2 for PG 1222+228 at a significance level of 
3.3cr (Mushotzky 1999). If we assume a "normal" power-law slope of a = —1.6 over the 0.2 to 2 
keV ROSAT band (Brandt et al. 1999), we find F x = 1.5 x 10~ 31 erg s" 1 cm~ 2 FLr 1 at 2 keV 
rest energy. Optical spectrophotometry (Wampler & Ponz 1985; Bechtold et al. 1984) implies 
F Q 1.1 x 10~ 26 erg s _1 cm~ 2 Hz -1 at rest wavelength 3000 A. (These are observed fluxes at the 
wavelength corresponding to the indicated rest wavelength.) From this, we find otox — -1.8 for PG 
1222+228. This result is uncertain because of the marginal X-ray detection and the possibility of 
variability, but at face value it is more consistent with a non-BAL than a BAL QSO. (Wilkes et 
al. 1994 quote an uncertain value 1.8 from Einstein data.) 

We conclude that the flux drop at A750 in PG 1222+228 is unlikely to be a BAL feature. 

2.2. An Instrinsic Lyman Limit System? 

Impey et al. (1995) suggested that the flux drop at A750 was an intervening LLS. We show 
below that Lyman limit absorption does indeed give a good fit to the spectrum. This fit, however, 
leaves open the question of the location of the absorbing gas. Because of the coincidence with the 
polarization rise, we consider here the possibility that the feature is an intrinsic LLS, associated 
with a high velocity outflow from the QSO. Impey et al. (1996) identify La absorption lines at z 
= 1.4857, 1.5238, 1.5272, and 1.5650 that might be associated with the A750 feature, if it is taken 
to be a LLS. A redshift of 1.486 corresponds to an outflow velocity ~ 60,000 km s -1 . As noted 
above, this is not unprecedented for a QSO outflow producing absorption lines. However, the lines 
associated with the redshift systems in question in PG 1222+228 are narrow, and narrow lines with 
relative velocities greater than 5000 km s _1 usually are assumed to be intervening. For a Lyman 
edge optical depth of unity, the measured equivalent widths of the La lines are consistent with a 
Doppler parameter b <J 30 km s _1 (see below), normal for an intervening LLS. In contrast, the 
mini-BALs observed at such high velocities have widths of order 1000 km s _1 (Hamann et al. 
1997). Could the feature in PG 1222+228 nevertheless be caused by ejected material with a high 
outflow velocity and a small velocity dispersion? 
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We are not aware of any other case in which a high velocity LLS with narrow lines has been 
shown to be intrinsic. However, the existence of many intrinsic, narrow, high velocity absorption 
systems in QSOs has been proposed by Richards et al. (1999). These authors compare the incidence 
of absorption-line systems, per unit relative outflow velocity, for highly luminous QSOs with that for 
less luminous ones. In the velocity range 5000 to 75,000 km s _1 , they find that absorption systems 
have a substantially larger frequency in luminous systems. Since intervening systems should have no 
dependence on QSO luminosity (assuming that discovery systematics are accounted for), Richards 
et al. conclude that at least the excess number of systems in the high luminosity QSOs are intrinsic. 

A peculiar absorption line ratio in the z = 1.94 system in PG 1222+228 has been noted 
by Ganguly et al. (1998). These authors present high resolution spectra that show two narrow 
components, separated by ~ 17 km s _1 , with very different strengths of the Al II and Al III 
absorption lines. Photoionization models indicate that the component with strong aluminum lines 
must have an anomalously high abundance of aluminum. This is reminiscent of claims of unusual 
abundances in BAL QSOs, including an excess of aluminum (e.g., Turnshek et al. 1996; Junkkarinen 
et al. 1997; Shields 1997). The reality of these abundance anomalies is in doubt, because of the 
effects of partial covering of the continuum source (Arav 1997). However, the basic observation of 
anomalously strong Al lines may be a possible parallel between PG 1222+228 and the BAL QSOs, 
where outflowing gas is clearly present. If this is a hint that the narrow, z = 1.94 system may be 
intrinsic, perhaps it adds plausibility to the idea that the z = 1.486 system (or its neighbors) may 
also be intrinsic. 

What might be the geometry of an intrinsic LLS in PG 1222+228? In order to explain the 
polarization rise, the absorber would have to intercept the line of sight to the continuum source 
but not the scattering source. The latter is often attributed to a wind driven off the inner edge 
of a "dust torus" (Krolik & Begelman 1986). The location of this may be related to the dust 

1/2 

sublimation radius, ~ 0.2L 4 g pc, where L4Q is the bolometric luminosity of the central source in 
units 10 erg s _1 (Laor & Draine 1993). The absorbing material, at a velocity of ~ 0.2c, would be 
at a smaller radius in order not to cover the scattering source. Its high velocity suggests an origin 
at a small radius where the escape velocity is of order the observed outflow velocity. The escape 
velocity from a central mass of M is 0.2c at a radius 10 15 ' 6 Mg cm, where Mg = M/10 9 M . Let us 
assume the absorbing material is at a radius of <J 10 18 cm, large enough to obscure the ultraviolet 
emitting part of an accretion disk but not the scattering source. At the observed speed, the crossing 
time would be a few years or less. Thus, the material should change radius substantially in the 
seven years between the Sargent et al. (1988) and the Impey et al. (1996) observations, and one 
might expect some change of velocity. These authors, however, quote velocities for the z = 1.486 
and 1.524 systems that agree within Az ~ 0.001. This corresponds to a change in outflow velocity 
of less than ~ 100 km s _1 . Such precise stability seems difficult to achieve. (Note, however, the 
stability of narrow features within the BAL profiles of some QSOs [cf., Weymann 1997]). A further 
problem involves the narrowness of the absorption lines. The radius of the ultraviolet emitting part 
of the disk would be at least ~ 10 gravitational radii, or about 10 15 ' 1 Mg cm. If the absorber is at 
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a radius ^ 10 8 cm, then the line of sight through the absorber to different parts of the continuum 
source would likely give noticeably different projected flow velocities. The observed linewidths (see 
below) are <; 30 km s _1 , less than one thousandth of the outflow velocity. 

These difficulties with the intrinsic absorber model for PG 1222+228 encourage us to examine 
the straightforward idea of an intervening LLS. 

3. INTERVENING ABSORPTION 

Impey et al. (1995, 1996) attributed the flux drop at observed wavelength A2300 to a LLS 
associated with the z = 1.4857 absorption line system. They identify absorption lines of H I, C III, 
N I, N II, Si II, and Si III. They also identify systems at 1.5238, 1.5272, and 1.5650, which have 
multiple Lyman lines and, for z = 1.5238, C III. Sargent, Steidel, & Boksenberg (1988) measure C 
IV AA1548, 1551 in the 1.486 and 1.524 systems with equivalent widths ~ 0.7 A. However, Steidel 
& Sargent (1992) give a spectrum showing no detectable Mg II absorption. 

Simple estimates suggested that the A2300 flux drop might be too gradual to be attributed to 
the converging Lyman lines of the z = 1.486 system. Therefore, we computed a model spectrum 
involving an assumed power-law continuum and absorption by the hydrogen lines and bound-free 
continuum. The column density of H I was parameterized by the Lyman edge optical depth, th, 
and the lines were assumed to have a Gaussian profile with a Doppler parameter set to a typical 
value b = 30 km s _1 . (A larger line width would give an excessive equivalent width for La for 
the required th • The observed decrement of the Lyman line equivalent widths actually suggests 
a narrower core line width, with some of the La equivalent width being attributable to a broader 
component of modest column density. These details do not affect our conclusions.) Inclusion of 50 
Lyman lines proved more than adequate to trace the convergence to the continuum optical depth. 
The model spectrum was convolved with a single component Gaussian instrumental line profile 
following the discussion of Impey et al. (1996), using a FWHM of 2.9 A for G190H and 4.2 A for 
G270H. If only the z = 1.486 system contributes to the LLS, then the Lyman line convergence is 
too close to the Lyman limit (A912) to fit the observed feature. However, allocation of some H I 
column density to both the 1.486 and 1.524 systems gave a good fit. Two additional LLS appear to 
be associated with the z = 1.938 and 1.174 systems. Figure 3 shows the resulting model spectrum, 
along with the observed flux, binned in intervals of ~ 7.3 A. The model has th = (1.0,0.5,0.8,0.6) 
for z = (1.174, 1.486, 1.524, 1.938), respectively. (The value of th for z = 1.174 is uncertain because 
of the poorly determined amount of scattered light below 2000 A.) The three distinguishable LLS 
show good agreement with the expected behavior of the Lyman continuum optical depth above 
threshold, for an intrinsic continuum slope in the range a ~ —1.5 to -2.0. This is consistent with the 
slope -1.8 found by Zheng et al. (1997) for their composite QSO spectrum in the wavlength range 
A600 to A1050. A value a = —1.8 is assumed in the fit shown in Figure 3. The gradual descent of 
the observed flux toward the Lyman limit for the z = 1.938 system is a puzzle, but it may involve 
the effects of unrelated absorption lines. An understanding of this is important, as it would play a 
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role in the classification of PG 1222+228 as a candidate Lyman edge QSO. Observations at higher 
spectral resolution would help to clarify the situation. 

Sargent, Steidel, and Boksenberg (1989) discuss the statistics of LLS in QSOs. For the redshift 
range in question, they give a mean incidence of LLS of N(z) ~ 1.5 per unit redshift. Their 
data show many QSOs with multiple LLS, although the number of LLS in PG 1222+228 may be 
somewhat higher than typical. However, the efforts to measure Lyman continuum polarization in 
QSOs to some extent targeted the candidate Lyman edge QSOs, and objects with LLS of moderate 
optical depth may have an enhanced probability to be included. 

We conclude that an intrinsic power-law continuum, together with cosmologically intervening 
Lyman limit absorption, provides a straightforward explanation of the ultraviolet spectrum of PG 
1222+228. 

4. THE INTRINSIC POLARIZED CONTINUUM 

The various LLS in PG 1222+228 substantially attentuate the observed continuum. What is 
the behavior of the polarized flux when corrected for the absorption? In view of the uncertainties in 
the measured polarization, a sufficient procedure is to estimate the polarized flux by multiplying the 
measured polarization in the chosen wavelength bins by the assumed intrinsic continuum flux. For 
this, we use the same wavelength bins described earlier, and the I V PL oc is^ 1 - 8 power-law continuum 
used in our fit. The resulting rotated Stokes flux, = q' x I\ PL , is shown in Figure 4. We see that 
the Stokes flux now rises strongly with decreasing wavelength in the region of the polarization rise. 
This resembles the result found for PG 1630+377 by Koratkar et al. (1995). A rising polarized 
flux is an important constraint on models for the origin of the polarization rise. 

SWH showed that the wavelength dependence of the flux and polarization in PG 1222+228 
and PG 1630+377 could be fit with an ad hoc model involving an accretion disk. The disk radiates 
as a black body, but the brightness temperature is depressed below the effective temperature 
for wavelengths below the Lyman limit, simulating a Lyman edge in the disk atmosphere. The 
polarization is assumed to rise abruptly at A912 by an arbitrary amount. Relativistic effects give 
a blueshifted, but still fairly abrupt polarization rise in the observed spectrum. The models are 
characterized by a*, the dimensionless angular momentum of the hole; the black hole mass; the 
accretion rate, Mq = M/(l Mq yr^ 1 ); and the viewing angle, iJ, b s = cos(# fc s ). SWH found that 
a* = 0.5 gave approximately the observed wavelength for the polarization rise, for a relatively 
edge-on viewing angle. Their fits to both objects had a fairly low value of T max , the maximum 
disk effective temperature, as required by the dropping flux in the Lyman continuum region. The 
correction for LLS absorption in PG 1222+228 hardens the far ultraviolet spectrum, and a higher 
value of T max is required to fit the energy distribution. Figure 5 shows the continuum flux and 
polarization for a model with a* = 0.5, Mg = 8.8 and Mq = 86. (We have assumed Ho= 70 
km s _1 Mpc _1 and qo= 0.5.) The model agrees reasonably well with the corrected flux in the Lyman 
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continuum and with the longer wavelength measurements. Although the corrected flux was assumed 
to be a power law, a disk continuum would also likely fit the observed flux, given some freedom 
to adjust the LLS optical depths. This model has a step-function rise in polarization from negligle 
polarization at wavelengths longward of A912 to an ad hoc value of 2.1 times the Chandrasekhar 
(1960) value for a pure scattering atmosphere. The model predicts an observed polarization and 
polarized flux that rise at a wavelength substantially blueshifted from A912, but the rise is more 
gradual than observed. This underscores the need for improved polarization measurements of 
this object. If the relativistic transfer function gives too gradual a polarization rise, even for an 
instantaneous polarization rise in the rest frame of the orbiting gas, then accretion disk models for 
the polarization rise will face a serious problem. 

The adopted model parameters give a bolometric luminosity L/Le = 0.36, where Le is the 
Eddington limit. Such a high value of L/Le is barely consistent with a geometrically thin disk. A 
larger value of a* would allow a larger mass for the required T max , but then the polarization rise 
would occur at a wavelength shorter than observed (cf. SWH). Evidently, an accretion disk fit to 
the corrected continuum of PG 1222+228, in the manner of SWH, pushes the disk parameters to 
the limits. Conceivably, the thickening of the disk corresponding to the large value of L/Le may 
be related to the origin of the polarization rise. 

5. DISCUSSION 

The Lyman continuum polarization rises are among the more puzzling recent observational 
discoveries concerning QSOs. The wavelength dependence, rising rather abruptly at nearly the 
same rest wavelength in the several known cases, suggests a connection with the bound levels of 
atoms. The proximity of the feature to A912 further suggests an association with the Lyman edge 
of hydrogen. The ad hoc model of SWH supports an association with the Lyman edge and raises 
the possiblity of confirming the presence of a relativistic disk and constraining its parameters. 
However, attempts to fit the feature with a physical model have encountered difficulties. This is 
an important problem for QSO theory. 

Are the reported polarization rises real? The coincidence of the polarization rise in PG 
1222+228 with a sharp drop in flux might raise the question of background problems with the 
FOS spectropolarimeter. Impey et al. (1995) argue that the polarization is unlikely to be less 
than 2.7% around 2000 A under any reasonable assumption for the FOS background. However, 
the degree of polarization in the far ultraviolet is uncertain by at least a factor two because of 
systematic errors involving background and scattered light in the FOS. The observed polarization 
rises in several QSOs occur at different observed wavelengths but similar rest wavelengths. We are 
not aware of any polarization rises of this nature in FOS spectropolarimetry of BL Lac objects 
or stars. There is an urgent need for a renewed capability for ultraviolet spectropolarimetry from 
space to confirm and extend the measurements. 
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Lyman continuum polarization rises have heretofore been associated with the candidate Lyman 
edge QSOs (see discussion by Koratkar et al. 1998). Our results suggest that PG 1222+228 
may not be a true member of this class. Measurements of the Lyman continuum polarization in 
additional QSOs are needed to clarify the frequency of occurence of the phenomenon and to look 
for correlations with features in the continuum flux, in the line intensities and polarization, and 
other properties. Observations to shorter rest wavelengths are needed to determine whether the 
polarization falls or continues to rise. Measurements of the time dependence of the polarization 
rises would be most interesting. The emitting radius of an accretion disk would be light weeks. 

Lee and Blandford (1995) considered a model for the far ultraviolet polarization rise of QSOs 
that did not involve the Lyman edge. Noting that a number of resonance lines of heavy elements fall 
in the rest wavelength range where the polarization rises, they suggested that resonance scattering of 
the QSO continuum might produce the observed polarization. Such a model could produce a rising 
polarized flux, since the polarization could be essentially zero at wavelengths without scattering 
contributions. As noted above, the polarization rise in PG 1222+228 may be too steep for models 
involving an accretion disk. In this case, alternative models such as resonance scattering may hold 
promise. We note that the polarized flux spectrum of PG 1630+377 (Koratkar et al. 1995) shows 
a strong rise at the wavelength of the N V emission line. 

The claim by Richards et al. (1999) that luminous QSOs have many intrinsic, narrow, high 
velocity C IV absorption systems has important implications. This would complicate the use of 
such systems to probe the evolution of galaxies and the intergalactic medium. The mechanism 
for producing the absorbing clouds would add another challenge to the subject of outflows from 
QSOs. Surveys of QSOs at different luminosities, to a uniform standard of signal-to-noise ratio, 
would allow confirmation of the claimed higher incidence of absorption in more luminous QSOs. 
Tests of the intrinsic nature of narrow absorptions have been summarized by Barlow, Hamann, &; 
Sargent (1997). These include time variability of the depth and profile of the lines and evidence 
for saturated but nonblack line profiles. Detection of changes in the velocities of the lines would 
be most revealing. Chemical abundances may be an indicator, given evidence for high abundances 
in the broad absorption and emission-line regions of QSOs (e.g., Hamann & Ferland 1993). If C 
IV systems in high and low luminosity QSOs have similar, mostly subsolar abundances, this might 
argue for their intervening nature. 

The author is grateful to C. Impey and C. Petri for providing advice and the reduced HST 
observations and to C. Sneden for the use of a computer subroutine. The work has benefitted 
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upon work supported by the Space Telescope Science Institute under Grant No. GO-07359.02. 
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Captions for Figures 

Fig. 1.— Ultraviolet spectrum of PG 1222+228 observed with HST by Impey et al. (1995, 1996). 
Figure gives flux at earth in 1CT 15 erg s^ 1 cm -2 A -1 as a function of observed wavelength. Vertical 
lines give the positions of the Lyman limit (above) and La (below) for redshifts z = 1.174, 1.486, 
1.524, 1.527, 1.565, and 1.938. Also shown is the Lyman limit for the emission-line redshift of 2.046. 
Data plotted here are from the G190H spectrum for A < 2224 and G270H for longer wavelengths 
(see text). Data is binned by 2 pixels. Data courtesy of C. Impey and C. Petri (1999). 

Fig. 2.— Flux and polarization of PG 1222+228 from HST observations by Impey et al. (1995, 
1996), binned as described in the text. Abscissa is rest wavelength in terms of the emission-line 
redshift. / is the measured flux in units of 10 -16 erg s _1 cm~ 2 A -1 . Rotated Stokes flux Q, called 
Q' in the text, is referred to the mean position angle of 168 degrees. The fractional rotated Stokes 
parameter, q = Q/I, is given as a percentage. 

Fig. 3. — Observed flux of PG 1222+228 compared with model involving power-law continuum 
and absorption by H I Lyman lines and continuum at four redshifts as described in the text. 
Observations have been binned into intervals of approximately 7.3 A for clarity. Axes are same as 
in Fig. 1 

Fig. 4. — Flux and rotated Stokes flux of PG 1222+228 before (filled squares) and after (open 
squares) correction for absorption by Lyman limit systems as described in the text. Also shown 
(crosses) are observations of PG 1630+377 by Koratkar et al. (1995). For PG 1630+377, the 
Stokes flux has been rotated to the mean position angle of 127 degrees, following Koratkar et al. 
The plotted values of I and Q have been scaled to correspond to the peak values in PG 1222+228 
in order to emphasize the wavelength dependences. The corrected, polarized flux of PG 1222+228 
(Q* in the text) rises strongly below A750, resembling the case of PG 1630+377. Axes are same as 
in Fig. 2. 

Fig. 5. — Accretion disk model compared with observations of PG 1222+228 (see text). Observa- 
tions are from Impey et al. 1995 (squares), Bechtold et al. 1984 (open circles), and Webb et al. 
1993 (open triangle). Model disk with a* = 0.5, M$ = 8.8 and Mq = 86 is viewed at angle [i b s = 
0.25. Axes are same as in Fig. 2. 
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Table 1. Continuum Polarization of PG 1222+228" 



Grating 




^max 


I 


?(%) 


u{%) 


P(%) 


e 


q' 


G190H 


1994 


2224 


4.04+0.04 


5.8±1.8 


-3.4+2.6 


6.4+2.0 


165±16 


6.7±1.9 


Both 


2224 


2287 


2.94+0.04 


2.4+2.0 


-3.8+2.1 


4.0+2.0 


151±26 


3.7+2.0 


Both 


2287 


2319 


4.59±0.05 


0.2±1.8 


-3.4+1.9 


2.9+1.9 


136±31 


1.5±1.8 


G270H 


2319 


2492 


9.27+0.04 


1.0±0.6 


-1.0+0.6 


1.2±0.6 


157±26 


1.3+0.6 


G270H 


2492 


2761 


9.33+0.02 


0.5+0.4 


-0.0+0.4 


0.3+0.4 


179+48 


0.5+0.4 


G270H 


2761 


3029 


14.14±0.03 


0.9+0.4 


-0.8+0.4 


1.2+0.4 


159±17 


1.2+0.4 


G270H 


3029 


3295 


16.85±0.04 


1.7±0.4 


0.2+0.4 


1.6+0.4 


4+14 


1.5+0.4 



a Data of Impey et al. (1995, 1996) rebinned with respect to the apparent LLS at A2300. Wavelength 
limits (observed) are given in A. Stokes parameters are given as I, the measured flux (units 10~ 16 
erg cm" 2 s _1 A^ 1 ), q = Q/I, and u = U/I. The quantity q' is Q/I rotated to P. A. 168°. Quoted 
uncertainties are la, derived from uncertainties in individual pixel measurements. Polarization p is 
corrected for bias (Wardle & Kronberg 1974, eq. A3). 



- 16 - 



20 




15 




10 




5 







-i 



i i i i I i i i i I i i i i I i i i i I i i i 



PG 1222+228 



J I I I I I I I I I I I I I I I I I I I I I I I I L 



700 800 900 1000 1100 1200 




Of 



4 
3 
2 
1 




~\ I i i i i I i i r 



J I L 



J I L 



700 800 



900 1000 1100 1200 
\ est (A) 



-17- 




-18- 



20 
15 
10 
5 




-i i | i i i i | i i i i | i i i i | i i i i | i i i n 

- D n □ □ □ x x x x * 

! x x x x >< x x x x x - 

r ■ x ■ 

: X X x X X 

■ 

1] I I I I I I I I I I I I I I I I I I I I I I I I I L 



700 800 900 1000 1100 1200 



20 
10 



i i i i i i i i i i i r 



i i i i | i i i r 

■ 1222obs 
x 1630obs 







fx ■ 



X | A X X * X x 



X - X X 



A X X 



d I I I I I I I I I I I I I I I I I I I I I I I I I L 



700 800 900 1000 1100 1200 



10 
5 



□ 1222corr 



~- 



.. , x X X x x 

J I I I I I I I I 



-x x— x- 



J I L 



X 



700 800 900 1000 1100 1200 

\ est (A) 



